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Abstract 


The optical radiation from the plasma discharge 
of an electron bombardment mercury ion thruster was 
investigated. This work extends and refines earlier 
measurements , correcting certain ambiguities that 
arose in the analysis of line amplitudes. Using 
the measured ratio of the Hg I line amplitude at 
3655 A to that at 3650 X, a theory incorporating a 
bimodal electron distribution (Maxwell electrons 
plus primary electrons) was used to obtain the av- 
erage electron temperature and primary electron 
fraction in the thruster ion chamber. The electron 
temperature ranged from about 1.2 eV to 6.6 eV ; 
whereas the primary electron fraction varied from 
zero percent to about 5 percent . These values de- 
pended upon the discharge voltage and the radial 
location of the measurement. The percentage of 
doubly ionized mercury produced in the chamber was 
also determined as a function of discharge voltage. 

Introduction 


ference and a low temperature Maxwellian distribu- 
tion. Although an attempt is made in Refs. 6 and 7 
to account for this bimodal distribution , there 
exists no satisfactory theoretical explanation to 
date. In order to interpret the observed spectral 
line amplitudes in terms of the electron energy 
producing the emitted radiation, it is convenient to 
incorporate this bimodal electron distribution in a 
theoretical model of the discharge plasma. Because, 
in principle, both primary and Maxwellian electrons 
can contribute to excitation and ionization of 
mercury atoms, a theoretical model of the discharge 
plasma should incorporate excitation and ionization 
by both groups. 

Theory 

The rate of populating the excited state 

due to collisional excitation from the ground state 
is given by 


In a previous paper (!) the authors presented 
the results of preliminary measurements of the op- 
tical spectra emanating from thruster plasmas. The 
present work extends and refines these measurements , 
correcting certain ambiguities arising in the earl- 
ier analysis. Basically, a method is described 
which permits evaluation of thruster discharge 
parameters from measured line amplitudes . 

The observed optical spectral amplitude is a 
measure of the number of emitted photons of energy 
corresponding to the wavelength of the spectral 
line. The problem is to relate this amplitude (pro- 
portional to the number of transitions at a given 
wavelength) to the physical processes populating 
the energy level producing the photon emission(2). 
The following statements can be made about the phy- 
sical nature of the discharge plasma within electron 
bombardment thrusters . 

The discharge operates at low pressure (less 
than 10-3 torr) so that collisional broadening of 
spectral lines, collisional energy transfer between 
atoms, and three body volume recombination would be 
expected to have a negligible effect on the observed 
line amplitudes (3) . With the possible exception of 
the resonant states, absorption should play a negli- 
gible role in excited state populations. 

The ion chamber plasma, although essentially a 
steady-state discharge, was not in thermal equili- 
brium at the electron temperature. This can be 
shown by comparing calculated electron temperatures 
for an equilibrium plasma(^>5) with Langmuir probe 
measured temperatures^, 7, 8>9) 
peratures 

higher than those calculated from equilibrium 
theory. The discrepancy results from application 
of "thermal equilibrium equations to nonequilibrium 
sustained plasmas(l°). 


The measured tem- 


Langmuir probe measurements of the energy dis- 
tribution of such discharges ' using a derivative 
method first described by Medicus(3l) indicate 
that the electron distribution consists of nearly 
monoenergetic , or primary, electrons at an energy 
detei'mined by the discharge chamber potential dif- 


= N o [N m S j (T e ) + Vj ( V lcm " 3 Se ° _1 (1) 

where Sj(T e ) is the Maxwell averaged excitation 
coefficient for the jth state and Sj(Ep) is the 
excitation coefficient averaged over the primary 
electron distribution (represented as a delta- 
function at energy Ep). 

The number of j -*• k transitions is then given 
by the steady-state relation 

M jk = ^j A jk cm ^ sec-1 (2) 

where Aj^ is the relative transition probability 
for the j ■* k transition, Ajj^ = Aj^/^Aj^* Equa- 
tion (2) can be rewritten with the aid of Equa- 
tion (l). 


N., = N N Cl + <f»X (T ,E )]S„ (T ) 
Jk o m Jk e’ p Jk e 
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where 
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and Sjfc = <QjkV e >, the brackets denoting averaging 
over the appropriate electron distribution. The 
quantities are related to the excitation cross 

sections Qj by the relation 


jk 
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They are the optical excitation functions for pro- 
ducing the observed spectral line corresponding to 
the j -*■ k transition. The most recent measure- 
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ments of optical excitation functions for mercury 
are those of Anderson, et al.( 12 ). A complete com- 
pilation of available mercury atom optical excita- 
tion functions can be found in Ref. 13. The num- 
ber of transitions given by Equation (2) is propor- 
tional to the spectral line amplitude N\ meas- 
ured by the photodetector of the monochromater . 

In order to obtain the electron temperature 
and primary electron fraction, it was necessary to 
ratio two line amplitudes. Examination of the ob- 
served Hg I mercury spectrum resulted in the selec- 
tion of spectral lines at 3650 X (63D2-63p 2 ) and 
3655 X (63 d2-63p 2 ). Population of the 63 d 3 level 
from higher lying states by cascading is negligible 
in the electron energy range of interest* 12 ^ and 
there is no singlet mixing of this state. The 
63Dp, however, mixes strongly with the 6 1 Dg to 
produce an optical excitation function that is 
characteristic of the singlet state^ 12 ). The ex- 
citation cross sections for these two lines are 
thus sufficiently different in their dependence on 
electron energy that the line ratios vary signifi- 
cantly. Thus a method based on that described in 
Refs. 10 and lU was used to obtain electron tem- 
perature and primary electron fraction. 

From Equation (3) the ratio of the number of 
3655 X transitions to 3650 X transitions is given by 

M 3655 = 

“3650 


1 + * X 3655 ( V E p } S 3655 (T e ) 

. 1 + * X 365 0 (T e’V WV. 


The variations of the appropriate excitation coef- 
ficients with electron temperature or primary elec- 
tron energy are given in Figs. 1 and 2. 

In order to use Equation (7) to obtain primary 
electron fractions <(> and electron temperatures 
T e , it was necessary to choose a value for the pri- 
mary electron energy Ep. Previous studies of the 
discharge of such thrusters (15, 16) indicated that 
the primary electron energy can be related to the 
discharge chamber potential difference by the ex- 
pression . 


E p = AVj - V p 


(8) 


where Vp is the plasma potential in the cathode- 
magnetic pole piece region. In the present study 
this potential, Vp, was taken to be 16 volts, based 
on electric probe measurements reported in Ref. 16. 
As noted in this reference, however, Vp is depend- 
ent on cathode geometry, and thus the value of 
16 volts is not universal. 

A qualitative estimate of the sensitivity of 
the transition number ratio ^36515/^3650 to Pri- 
mary electron energy (and thus the value of Vp 
used in Eq. (8)) is shown in Fig. 3. Here the left 
side of Equation (7) was plotted against Ep at 
two electron temperatures and three values of pri- 
mary fraction. It was evident from such plots 
that the transition number ratio became less sensi- 
tive to changes in primary energy as the primary 
fraction decreased and the electron temperature in- 
creased. Thus, an error in the selection of Ep 
would have the largest effect at the lower electron 
temperatures and high primary fraction. 

In summary, the theoretical model developed 
herein incorporates excitation by both Maxwellian 
and primary electrons. This represents a convenient 
model having simple features that allow for 
straightforward interpretation of the optical radia- 


tion from the thruster discharge chamber. Recent 
Langmuir probe measurements indicate that such a 
distribution does indeed exist in thruster dis- 
charges. The analysis uses measured spectral line 
amplitudes to calculate the electron temperature 
and primary electron fraction in the discharge cham- 
ber plasma of a bombardment thruster. Once these 
parameters are determined, other properties of the 
discharge such as ion fractions and density varia- 
tions can be studied. 

Experiment 

The preceding theory was used to interpret the 
results of spectroscopic measurements obtained on 
a 30-cm diameter hollow cathode bombardment thrus- 
ter. Detailed descriptions of such thrusters, 
shown schematically in Fig. 1*, can be found in 
Refs. 17 and lR. The hollow cathode serves as the 
source of ionizing electrons which are contained in 
the discharge chamber by means of field-shaping 
permanent magnets. Propellant utilization, un- 
corrected for doubly charged mercury ion production, 
ranged from about 80 percent to about 100 percent. 

The extraction grid open area was about 72 percent . 

An Elbert mounted, plane grating monochromator 
with a 0.5-meter focal length was used to study 
discrete spectral line amplitudes. The instrument 
dispersion was 16 X/mm. The plane grating could Q be 
rotated Jo cover the wavelength range from 2000 A 
to 8000 A. The recorded photomultiplier signal was 
corrected for the spectral response of the optical 
system. This response was determined using a ra- 
diance standard calibrated by the National Bureau 
of Standards . 

The thruster was operated in the same 7.6-meter 
diameter by 18. 3-meter facility described in Ref. 1. 
There was a question as to how much of an error 
would be introduced as a result of viewing the ion- 
ization chamber through the exhaust beam. This 
downstream radiation was measured by locating the 
monochromator normal to the thruster axis and view- 
ing a region Just downstream of the thruster. Such 
measurements indicated that the downstream radia- 
tion was about two to three orders of magnitude less 
than that produced in the discharge chamber and could 
thus be neglected. 

The monochromator was used to obtain radial pro- 
files of the 30-cm diameter thruster. The proce- 
dure consisted of setting the monochromator on the 
peak of a given spectral line. A motor driven 25 
micron horizontal slit was then used to scan the 
image on the vertical entrance slit of the mono- 
chromator. In this manner the area of thruster 
cross section observed represented about 0.005 per- 
cent of the total thruster cross sectional area. For 
measurements of the thruster discharge chamber 
plasma, profiles were obtained along the radial 
direction away from the neutralizer location as 
shown in Fig. 5. 


Results 


Spectral line amplitudes were obtained for dis- 
charge chamber potential differences AVt of 33, 

38, 1*3, and 58 volts. The discharge current was 
maintained at 8 amperes, and the thruster beam cur- 
rent was constant at 1.5 amperes. At each dis- Q 
charge vgltage the ratio of the number of 3655 A 
to 3650 A transitions was calculated at four ra- 
dial positions. These ratios are shown graphically 
in Fig. 6. Here the measured ratio l3655/™3650 is 
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plotted at the four radial positions . Of particu- 
lar interest is the minimum occurring at a primary- 
electron energy of 22 eV (AVj = 38 V) at all radial 
positions except that furthest from the thruster 
axis. It will he shown that such a minimum can 
he explained hy the existence of nonrandomized pri- 
mary electrons in the plasma discharge . 

In Fig. 7, Equation (7) is plotted against 
electron temperature for the four primary electron 
energies studied. Several curves corresponding to 
different values of primary electron fraction <ji 
ranging from 0 to 0.1 are presented. The problem 
is to correlate these curves with the experimental 
curves of Fig. 6. In this way, the 
ratio, measured as a function of primary elec- 
tron energy, can yield the electron temperature and 
primary electron fraction corresponding to a parti- 
cular primary energy. Note that in Fig. 7, 
curves of > 0 also show minima, analogous to 
the minima of the data curves of Fig. 6 at a pri- 
mary electron energy of 22 eV. It is this fact 
which suggests the procedure to be used to inter- 
pret the spectral amplitude data. In reducing the 
data, use is made of the following conditions: 

Condition (l): Increasing the primary electron 

energy increased the available discharge power per 
atom, presumably increasing the electron tempera- 
ture. Thus the electron temperature obtained from 
Fig. 7(d) would be greater than the temperature 
obtained from Fig. 7(a). 

Condition (2): Because an increase in primary 

electron energy increases the mean relaxation time 
for primary electrons (19) f it was assumed that the 
primary electron fraction did not decrease with in- 
creasing primary electron energy. The average pri- 
mary fraction obtained from Fig. 7(d) would thus 
be greater than the average fraction obtained from 
Fig. 7(a). 

Condition (3): For values of Ep up to 22 eV 

and for normalized radii less than 0.75, the meas- 
ured N3655/N3650 ratio (Fig. 6) decreased with 
increasing electron energy. Also, the theoretical 
N3655/N3650 ratios of Fig. 7 decreased with in- 
creasing electron temperature for temperatures up 
to about U eV. Thus for Ep 22 eV, portions of 
the ^-curves of Figs. 7(a) and (b) to the right 
of the minima locus (unused region) can be neglect- 
.ed. Similarly, for Ep > 22 eV, portions of the <(>- 
curves of Figs. 7(c) and (d) to the left of the 
minima locus (unused region) are neglected. 

The allowable range of values for primary 
electron fraction at different radial locations was 
obtained using Condition (2) and the following pro- 
cedure. For a given radial location (e.g., r = 0) 
and primary electron energy equal to 22 eV, a line 
was drawn parallel to the abscissa of Fig. 7(b) so 
as to intersect the ordinate at a value equal to 
the measured N 3655 /N 3650 ratio (at r = 0, the 
ratio was 0.065, from Fig. 6). This line is tan- 
gent to a particular <j>-curve at its minimum, and 
thus intersects the locus of minima at a particular 
value of <(> and T e . (For r = 0, the measured 
H 3655/N3650 line was tangent to the <P = 0.C3 
curve at a temperature T e of about 3-5 eV.) This 
value of 4> represents an upper bound of the <f> 

(and also T e ) range at the particular radial lo- 
cation and for primary electron energies from 17 eV 
to 22 eV. Lower bounds on 4> and T e for this 
primary electron energy range and for a particular 
radial location were obtained from the measured 
JI 3 g 55 /M 36 5 o ratio at 17 eV and represented by a 
horizontal line overlayed on Fig. 7(a), The inter- 


section of this line with the 4>-curves to the left 
of the minima and the assumption that for tenuous 
plasmas, the electron temperature is generally ex- 
pected to exceed about 1 eV'i served to obtain 
these lower bounds (Thus, for r = 0, where 
%655/^3650 = 0.074 at Ep = 17 eV, the lower bound 
on from Fig. 7(a) was found to be about 0.01 at 
an electron temperature of about 1.2 eV.) The 
lowest lower bound on <)> was obtained at a normal- 
ized radius of 0.75 and was <j> = 0. 

For Ep > 22 eV, the procedure used to obtain 
the appropriate ranges of electron temperature and 
primary electron fraction was as follows: The 

maximum primary electron fraction was obtained from 
the intersection of the horizontal line correspond- 
ing to the measured ^3655/^3650 with the. 

curves of Fig. 7(d). The largest value of this 
ratio (equal to 0.089) was measured at r = 0, for 
which the largest value of <p intersected was 
about 0.05. Thus, for the lower primary energy of 
28 eV (Fig. 7(c)), the intersection of measured 
N3655/N365O values with curves of <f> > 0.05 were 
disregarded in accordance with Condition (2). 

The minimum value of <|> at a given normalized 
radius and for Ep > 22 eV could not be less thah 
the value of 4> obtained at 22 eV, again in ac- 
cordance with Condition (2), Thus for r = 0 and 
Ep > 22 eV, the range of <j> was defined to be 0.03 
4> i 0.05. For Ep 22 eV and r = 0, the range 
was 0.01 <_ p 0.03. Using these 4> ranges and 
the appropriate curves of Fig. 7, the range of 
electron temperature for r = 0 and different 
values of Ep were also determined. Analogous pro- 
cedures to those Just described were used to obtain 
the ranges of 4> and T e at the other three radial 
locations studied. 

The results of the analysis are summarized in 
Figs. 8(a) through (d). Here the variation in the 
range of electron temperature with primary electron 
energy is presented for four radial positions. The 
boundaries of the darkened regions were deter- 
mined by the range of primary electron fraction. 

The values of electron temperature obtained were in 
general agreement with Langmuir probe .measurements 
of hollow cathode thruster discharges' ' . The 
fact that the largest values of d> and T e were 
obtained on the thruster axis, whereas the smallest 
values of these parameters were obtained at three- 
quarters radius, is consistent with recent meas- 
urements of a primary electron region centered about 
the thruster axis(20). 

Effects of metastables . - For the bombardment 
thruster discharge chamber plasma, the contribution 
to excited states by collisional excitation from 
metastable states can be shown to be much less than 
the direct excitation from the ground state and 
could thus be neglected. The two important meta- 
stable states in mercury are the 63P 0 and 63p 2 
states. The ratio between the rate of collisional 
excitation of the 6^3 state of mercury from the 
63p 0 metastable level to the rate of collisionally 
populating the 6^3 from the ground state was 
calculated to be about 0.005 over an electron tem- 
perature range of 1 to 7 eV. The corresponding 
ratio for excitation from the 63Pg metastable 
state was about 0.027 over the same electron tem- 
perature range. These ratios were lower at electron 
temperatures below 4 eV. The calculations were 
based on the Gryzinski theory( 21 > 22) and the as- 
sumption of a pure Maxwellian electron distribution. 
Experimentally , the primary electron density was 
generally less than 5 percent of the Maxwellian 
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electron density. This fact, in conjunction with 
the fact that excitation coefficients for primary 
and Maxwellian electron excitation are of the same 
order of magnitude (Figs. 1 and 2) implies that 
these estimates of the metastable level contri- 
bution to Equation (l) are reasonable. Thus, the 
exclusion of the metastable contribution to the 
collisional excitation of the 6^D states was jus- 
tified. 

Ion excitation . - For ion excitation the spec- 
tral amplitude at 3984 A, corresponding to the 
6 2 P5/2 - 6s 2 2 Dg/ 2 transition in Hg II, was used 
as the diagnostic line. This line produced the 
strongest signal in the near visible Hg XI excita- 
tion spectrum. The calculated excitation cross 
section for the resonant 6 2 P3/ 2 level is given in 
Fig. 9. From the work of Penkin, et al.( 23 ) one 
would expect cascading to contribute to the popula- 
tion of this level. Available data, however, are 
insufficient to estimate the extent of the cascade 
contribution. Also shown is the calculated level 
excitation cross sections for the 5d| /g 3/26s 2 (J = 
4) Hg III line at 4797 8. This line whs’ufeed to 
obtain the fraction of Hg III in the discharge 
plasma. These cross sections are plotted as func- 
tions of electron energy. Excitation from the 
ground state of the ion (single ion for the 3984 A 
line and double ion for the 4797 A line) was as- 
sumed. 'From comparisons between such calculated 
cross sections for helium and experimentally meas- 
ured helium optical excitation functions* 22 ' the 
calculated cross sections for mercury are expected 
to be within an order of magnitude of actual cross 
sections. 

Density variations . - The variation of the ion 
to neutral fraction was studied as a function of 
primary electron energy at different radial loca- 
tions in a 30-cm hollow cathode thruster. Using 
the method to obtain Equation (7) measured ampli- 
tudes at 3984 A for Hg II and 3650 A for Hg I were 
ratioed to obtain the expression 


5+ _ %984 1 + 4’X5650( E p a T e) s 565o( T e) 
N 0 H 3650 1 + ‘ t>X 3984^ E p’ T e^ S 3984^ T e^ 


where %9Q4(Te) 


( Q 6 2 R 


3/2 


r e) a 3984- 
/ m 


O) 


The quantity A 3 ggit is the relative transition 
probability . Its value is unknown, but is expected 
to be substantially less than one. This is because 
the 6 2 P ?) / 2 upper state of the 3984 A transition is 
a resonant state and thus would radiate to the 
ground state. The excitation coefficient ratios 
Xx(E p ,T e ) have been defined earlier. Because 
ffgggi, is constant, the right side of Equation (9) 
yields the relative variation of N+/No. The 
range of values of <t> and T e used were obtained 
from Fig. 8. The relative variation of N+/Nq is 
given in Table I (third column) as a function of 
Ep and radius . The calculations indicate that the 
ion fraction ’'4 /Nq did not vary appreciably with 
discharge voltage, showing a maximum at a primary 
electron energy of 22 eV for radial locations near 
the thruster axis. 

d'he relative variation of N + can be obtained 
from the relation 



where it is assumed that N+ a. N m - Combining Equa- 
tion (10) and Equation (9), the relative variation 
of Nq can also he obtained. These variations are 
shown in Fig. 10, in which an average of the normal- 
ized ion and neutral atom densities are plotted as 
functions of primary electron energy. Maximum 
densities occurred at the lowest energy (17 eV ) . 

For each radial location the data were normalized 
to unity at this energy. The marked decrease in 
both ion and atom densities can probably he at- 
tributed to the decrease in propellant flow re- 
quired to maintain the constant beam current 
(l.5 amperes) with increased discharge chamber po- 
tential difference. This fact suggests the pos- 
sibility of relating propellant flow changes to 
changes in atom and ion spectral amplitudes. Also 
noted was a leveling off of the neutral density 
at high discharge chamber voltages. This result 
is in agreement with a conclusion of Ref. 20, stat- 
ing that the neutral density remains constant at 
high discharge power levels (high eV/ion). 

Fraction of doubly charged ions . - Q Excitation 
of doubly-charged mercury ions at 4797 A were also 
measured in the 30-cm diameter thruster. An ex- 
pression analogous to Equation (9) (with U 79T A 
substituted for 398*1 X) can be used to obtain 
N++/Nq. In this case, = 1> because the upper 

state of this transition has only a single path for 
radiative decay* 2 ^'. Thus the spectral amplitude 
ratio gives the fraction of H++ in the discharge 
directly. These fractions are given in the last 
column of Table I. The accuracy of these measure- 
ments depends markedly on the accuracy of the ex- 
citation coefficients used. Because no measured 
excitation functions for Hg III are available to 
the authors ' knowledge , it is not possible at pre- 
sent to assess the accuracy of these calculations. 

It should he emphasized that the double ion frac- 
tions presented here cannot he directly compared 
with mass spectrometer measurements of double ions 
in the beam* 25 ). This is because the fraction 
N+/Nq is not known, and beam measurements yield 
N++/N + ratios . 

It should he noted that excitation of Hg III 
at an electron energy of 17 eV was observed at all 
radial positions, although quite weakly. The 
existence of Hg III at such low electron energies 
is to he expected because the threshold for ioniz- 
ing Hg II to Hg III from the 6 2 Dtj/g metas+.able 
level is about 15 eV . In addition, the "tail" of 
the low temperature Maxwellian distribution of 
electrons can produce both ionization of Hg II 
and excitation of Hg III. 

Concluding Remarks 

Spectral amplitude measurements indicated that 
the discharge was a nonequilibrium sustained plasma. 
The energy distribution of electrons in the dis- 
charge of hollow cathode thrusters can be repre- 
sented by a distribution of monoenergetic electrons 
superimposed on a Maxwellian electron distribution, 
as has been suggested by other workers. Spectral 
line amplitudes at 3655 A and 3650 K were ratioed 
to obtain electron temperatures and ratios of pri- 
mary to Maxwellian electron densities at four 
mary electron energies and different radial Loca- 
tions. It was found that the electron temperature 
ranged from about 1.2 to 6.6 eV ; whereas the primary 
to Maxwellian electron density ratio ranged from 
zero to about 0.05. These values deoended upor> pri- 
mary electron energy and radial location in the 
thruster. 
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Normalized ion and neutral densities were found 
to decrease with increasing discharge voltage at 
constant emission current and beam current. It 
was presumed that this behavior was attributable to 
a decreasing propellant flow with increasing volt- 
age. This decrease in flow was required in order to 
maintain constant beam current . The neutral density 
leveled off at high discharge voltages (high eV/ion) 
in agreement with earlier reported work. 

The fraction of doubly charged mercury ions 
did not exceed about 2 percent of the. atom density 
in the discharge plasma at a discharge voltage of 
58 V. At lower discharge voltages this fraction 
did not exceed about 0.8 percent. 

One of the primary purposes of this investiga- 
tion was to determine the extent to which the op- 
tical radiation emanating from the ion chamber of 
electron bombardment thrusters could be used to 
study the discharges of such thrusters. The high 
sensitivity of the radiation output to changes in 
thruster parameters suggests potential engineering 
applications. For example, studies are presently 
underway to investigate the light output from the 
neutralizer region. Such studies should lead to a 
better understanding of the operation of this thrus- 
ter component. It is possible that the light out- 
put could be used in a neutralizer controller. It 
is concluded that this study demonstrated the possi- 
bility for using spectroscopic techniques as diag- 
nostic tools for thruster research. 
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Figure 1. - Excitation coefficients for the b^-^Potransi- 
tion of Hgl at 3650 A. 
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Figure 4. - Section view of a 30-centimeter-diameter thruster. 
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Figure 5. - 30-Centimeter-diameter thruster image 
viewed from downstream of the extraction grids, 
72 percent open area. 



Figure 6. - Measured transition number ratio against 
primary electron energy, where E p * AVj - 16 eV. 
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Figure 7. - Theoretical transition number ratio against electron temperature. 
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Figure 10. - Normalized ion and atom density variations with elec- 
tron energy. 30-Centimeter-diameter hollow-cathode thruster; 
ion beam current, 1.5 amperes. 
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